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ABSTRACT 

We present 30 GHz measurements of the angular power spectrum of the cosmic microwave back- 
ground (CMB) obtained with the Sunyaev-Zel'dovich Array (SZA). The measurements are sensitive 
to arcminute angular scales, where secondary anisotropy from the Sunyaev-Zel'dovich effect (SZE) 
is expected to dominate. For a broad bin centered at multipole 4066 we find 60±f 5 ,uK 2 , of which 
26 ± 5 ijK 2 is the expected contribution from primary CMB anisotropy and 20 ± 28 fiK 2 is the ex- 
pected contribution from undetected radio sources. These results imply an upper limit of 149 /iK 2 
(95% C.L.) on the secondary contribution to the anisotropy, lower than that reported previously by 
other 30 GHz instruments. The SZA interferometric observations employed a hybrid configuration 
of antennas including short and long antenna separations to provide high sensitivity to arcminute 
anisotropy while simultaneously detecting 30 GHz radio sources at much higher resolution. The hy- 
brid configuration was also used to check whether SZE anisotropy power was suppressed by projecting 
out SZE power from clusters of galaxies correlated with radio sources. No significant effect was found. 
The level of SZE anisotropy power is consistent with expectations based on recent determinations of 
the normalization of the matter power spectrum, i.e., as ~ 0.8. 

Subject headings: cosmic microwave background - cosmology: observations - cosmological parameters 
- large-scale structure of universe - techniques: interferometric 



1. INTRODUCTION 

Density perturbations at the epoch of recombination 
are imprinted onto the Cosmic Microwave Background 
(CMB), leaving temperature anisotropy that has now 
been well-studied on a wide range of angular scales. On 
scales of several arcminutes and smaller, correspond- 
ing to multipole moments of £ > 3000, the level of 
CMB anisotropy power from primordial fluctuations is 
strongly suppressed by photon diffusion, and secondary 
sources of power, including the Sunyaev-Zel'dovich ef- 
fect (SZE), are expec ted to play a significant role (e.g., 
iHu fc Podelson1 l2002). The SZE results from the inverse 
Compton scattering of CMB pho tons by the hot elec- 
tron gas within clusters of galaxies (jSunvaev fc Zeldovichl 
1972). This interaction leaves a small spectral dis- 



of 



IL 



Electronic address: sharp@oddjob.uchicago.edu 

1 Kavli Institute for Cosmological Physics, University 
Chicago, Chicago, IL 60637 

2 Department of Physics, University of Chicago, Chicago, 
60637 

3 Department of Astronomy and Astrophysics, University of 
Chicago, Chicago, IL 60637 

4 National Radio Astronomy Observatory 

5 Enrico Fermi Institute, University of Chicago, Chicago, IL 
60637 

6 Owens Valley Radio Observatory, California Institute of Tech- 
nology, Big Pine, CA 93513 

7 Department of Space Science, VP62, NASA Marshall Space 
Flight Center, Huntsville, AL 35812 

8 Columbia Astrophysics Laboratory, Columbia University, New 
York, NY 10027 

9 Department of Physics, Columbia University, New York, NY 
10027 

10 Department of Astronomy, Columbia University, New York, 
NY 10027 



tortion in the CMB which produces anisotropy power 
on scales of I ~ 2000 — 10000, detectable as a decre- 
ment in the CMB intensity at 30 GHz. The amplitude 
of this power is extremely sensitive to the history of 
structure formation and, specifically, to the value of as, 
the normalization of th e matter power spectrum (e.g., 
iKomatsu fc Se liak 2002). Evidence for anisotropy above 
that expected from the primary CMB on these small an- 
gular scales has been detected at 30 GHz by the Cosmic 
Background Imager (CBI) and Berkele y Illinois Mary- 
land Association (BIMA) experiments (Rcadhca d et al.l 
2004!: iDawson et alJ feOOe). Recent constraints on ex- 
cess power at these scales from the Arcminute Cosmol- 
ogy Bolometric Array Receiver (ACBAR) at 150 GHz 
indicate that the spectrum of the reported excess is in- 
consistent with thermal CMB flu ctuations, but could b e 
consistent with SZE fluctuations (|Reichardt et al.| [2008). 
Taken together, these measurements indicate a level of 
SZE anisotropy power consistent with a value of as some- 
what greater than those preferre d by other contemporary 
measurements of the param eter (i Vocvod kin fc Vikhlininl 
[200l IKomatsu et al.|[200cl . 

In this paper, we describe a new, high-sensitivity mea- 
surement of power in the CMB on scales ranging from 
t ~ 3000 - 6000 made at 30 GHz with the Sunyaev- 
Zel'dovich Array (SZA). 

2. OBSERVATIONS 

2.1. SZA Data 

The CMB anisotropy data presented here were 
obtained with the Sunyaev-Zel'dovich Array (SZA), 
an eight-element interferometer located at Caltech's 
Owens Valley Radio Observatory near Bishop, Cali- 
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fornia. The SZA antennas are equipped with sensi- 
tive, wide-bandwidth receivers operating at 30 GHz 
and 90 GHz. For these observations we used the 
30 GHz receivers, tuned to detect sky frequencies of 
27—35 GHz. The receivers are based on low-noise, cryo- 
genic high electron mobility tr ansistor (HEMT) ampli- 
fiers (|Pospieszalski et ail 119951 ). with characteristic re- 
ceiver temperatures T rx ^ 11 — 20 K. Including atmo- 
spheric and other noise contributions, the typical system 
temperatures were 35—45 K. The 3.5-m SZA antennas 
have a primary beam that is well-described by a Gaus- 
sian with a FWHM of 12' at the center of the 30 GHz 
band. Cross-correlations of the signals from pairs of the 
eight antennas (visibilities) are formed in a digital cor- 
relator, which processes the 8 GHz IF bandwidth in 16 
bands of 500 MHz, each of which is further subdivided 
into 17 channels of 31 MHz, allowing rejection of narrow- 
band interference. 

For the anisotropy measurements reported here, the 
SZA antennas were arranged in a hybrid configuration 
to provide simultaneous sensitivity to the arcminute- 
scale structure of the SZE signal from galaxy clusters 
and the finer-scale contamination from r adio sources (de- 
tails of the configuration can be found in lMuchovei et all 
2007). Six of the eight antennas were packed closely to- 
gether (spacings of 4.5—11.5 m), yielding 15 baselines 
with typical projected lengths of 400— 1400A, correspond- 
ing to arcminute angular scales; the window function 
for our anisotropy measurement is determined by the 
baselines formed from combinations of these six anten- 
nas. By convolving the distribution of projected baseline 
lengths in the unfiagged data with the autocorrelation of 
the antenna illumination pattern, we obtain the ^-space 
window function shown in Figure [TJ This is the filter 
through which we observe the power spectrum; multi- 
plying this function by £(£ + 1)Ci/2tt and integrating 
over £ yields our measurement, where Ct is the CMB an- 
gular power spectrum. The mean value of the window 
is £ = 4066, with 68% of the area encompassed by the 
interval £ = 2929 — 5883. In addition to the central six 
antennas, two antennas were positioned ~50 m from the 
array center. The 13 baselines involving these antennas 
are sensitive to ~ 20" scales, corresponding to multipoles 
of £ ~ 30000, and can be used to identify radio sources 
that contaminate the short-baseline data. 

Between November 2005 and June 2007, we devoted 
1340 hours of observations to 44 SZA deep fields, com- 
prising 2 square degrees. These fields were selected to 
pass nearly overhead at the ~ 37° latitude of the tele- 
scopes, and arranged in right ascension to accommo- 
date other large SZA observing programs. Within these 
constraints, the fields were placed at locations with low 
Galactic emission in the IRAS 100 /an survey. No con- 
sideration was given to the presence of radio sources to 
avoid any bias that might be introduced by the correla- 
tion between radio sources and galaxy clusters. One of 
the 44 fields was found to contain a 700 mJy radio source 
(> 4000cr); unable to exclude this source from our data 
with sufficient precision, we are forced to omit this field 
from our analysis. Using simulated SZE sky maps we 
have verified that even in the pessimistic case in which 
this single field contains the most power of any of our 
fields, the bias introduced by excluding it is 10%. 

Observations were designed to allow subtraction of 
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Fig. 1. — (Upper Panel) CMB anisotropy power measured at 
30 GHz by the SZA (blue circle), CBI (red triangles), and BIMA 
(green square) experiments. Also shown are the WMAP5 primary 
power spectrum (dashed) and a model 30 GHz SZE spectrum (dots) 
for CT8=l-0 (upper) and 0.8 (lower). The sums of the primary and 
secondary power spectra are shown as solid lines. All data are 
corrected for contributions from undetected radio sources. (Lower 
Panel) For the same experiments, we show the window functions 
of the measurements, normalized to their peak values. 

ground contamination. Fields were observed at constant 
declination in groups of four, spaced by four minutes of 
right ascension. We observed the fields in 16-minute 
blocks between observations of a nearby phase calibra- 
tor. Each field, starting with the westernmost field of 
the group, was observed for just under four minutes at a 
time, so that after including time lost to slews and cal- 
ibration we tracked all four fields through precisely the 
same azimuth/elevation path. This ensures that each 
sees the same contribution from the ground, which can 
be removed by subtracting the mean of the four fields 
(in practice we implement this subtraction by means of 
a constraint matrix, described in §[3]). At an integration 
time of 20 seconds, we obtain 10 integrations per field in 
each cycle. The total integration time per field was ~ 25 
hours. Of this, approximately 20% of the data were dis- 
carded due to hardware problems, antenna shadowing, 
excessive noise, or jumps in calibrator phase or system 
temperature which occurred on a time-scale faster than 
our calibrations. An additional ~5% of the data were 
flagged for showing unexpected correlations among base- 
lines and bands or statistically unlikely noise behavior 
on a single antenna. These cuts significantly improved 
the results of the jackknife tests discussed in § 13.3) the 
data quality cuts were refined using only the jackknifed 
datasets to avoid biasing ourselves against signal. The 
first member of each group of four SZA fields is listed 
in Table [TJ along with its position, the total, unfiagged, 
on-source integration time, and the achieved RMS noise 
level for the field, both with and without data from the 
long baselines. 

The absolute calibration of the flux density scale is de- 
rived frombimonthly observations of Mars. We use the 
iRudvl (|1987t ) model to predict the brightness tempera- 
ture of Mars as a function of frequency and time. Since 
the planet is partially resolved by our longest baselines, 
a strong, unresolved source is used to transfer the cali- 
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TABLE 1 
SZA Anisotropy Fields 



Group of 


1 st Field RA 


Declination 


integration time 


RMS noise (short baselines) 


RMS noise (all baselines) 


4 Fields 






[hrs] 


[mjy] 


[mjy] 


cmb07 


02:07:37.00 


+34:00:00.0 


22.9 


0.18 


0.12 


cmbAl 


02:12:00.00 


+33:00:00.0 


19.7 


0.21 


0.14 


cmbll 


02:12:00.00 


+32:37:08.2 


22.8 


0.21 


0.14 


cmbRl 


02:12:15.60 


+32:11:24.8 


32.3 


0.18 


0.12 


cmbYl 


02:12:00.00 


+31:51:24.4 


23.0 


0.18 


0.12 


cmbEEl 


14:18:39.24 


+35:31:52.3 


20.9 


0.19 


0.14 


cmbXXb 


21:28:50.60 


+24:59:35.0 


18.2 


0.21 


0.15 


cmbDDl 


14:18:40.10 


+35:01:42.0 


18.0 


0.26 


0.17 


cmbAAl 


21:24:38.70 


+25:29:37.0 


19.6 


0.22 


0.16 


cmbCCl 


02:11:31.30 


+33:27:43.0 


21.9 


0.19 


0.14 


cmbBBl 


21:24:38.10 


+25:59:24.0 


18.8 


0.21 


0.15 



bration from the short baselines. The absolute calibra- 
tion was cross-checked by comparing SZA observations 
of Jupiter to those of W MAP and CBI (jPage et al.ll200l 
iReadhead et alj |2004[ ) . Based on these measurements, 
we conservatively estimate that our absolute flux scale is 
accurate to better than 10% (20% in power). 

2.2. VLA Data 

Away from the Galactic plane, the most significant 
contributor of non-CMB power on arcminute angular 
scales at 30 GHz is compact radio sources. Since the 
power from these sources is constant as a function of £, 
their contribution to l[l + l)Ci/2n can be quite large 
at small angular scales. While the brightest sources at 
30 GHz can be detected by the long baselines of the SZA, 
sources near our noise level (and hence undetectable in 
our data) can still contribute substantial power to the 
measurement. We therefore supplement our data with 
a higher-sensitivity search for radio sources in our fields 
using the NRAO's Very Large Array (VLA). 

The large disparity between the VLA (25 m) and SZA 
(3.5 m) antenna sizes makes it impractical to survey the 
SZA fields at a frequency near 30 GHz. Instead, we 
elected to conduct the survey at 8 GHz, a compromise 
between survey speed, which decreases toward higher fre- 
quency as the VLA primary beam shrinks, and spectral 
extrapolation of the detected sources. In November 2006, 
we collected 36 hours of data with the VLA in its C con- 
figuration, targeting the 24 fields that the SZA had ob- 
served at that point. An additional set of observations, 
this time in the more compact D configuration, which 
better matches the angular resolution of the SZA long 
baselines, covered all but two of the remaining fields in 
August 2008. Because the VLA primary beam at 8 GHz 
is only 4.9' across, it was necessary to make mosaic obser- 
vations with 19 pointings to cover the area of the SZA's 
larger primary beam. 

We reached a noise level of ~ 30 ^Jy at the center 
of this mosaic pattern, less than 1/6 of the noise level 
in our SZA data. The integration time is tapered at 
pointings away from the center of the mosaic, such that 
the achieved sensitivity in the 8 GHz maps has a profile 
that matches that of the SZA at 30 GHz, determined by 
our primary beam pattern. 

The VLA data were calibrated and mosaicked in AIPS. 
The map sensitivity a was determined from the mosaic 
images using the task RMSD. Sources were then ex- 
tracted from the maps using the task SAD, with a lim- 



iting significance of 5er, yielding more than 180 sources. 

Due to weather and scheduling limitations, two of the 
SZA fields were not observed with the VLA at 8 GHz. 
These fields fall within the cover age area of the 1.4 G Hz 
NRAO VLA Sky Survey fNVSS [Condon et al.lll99l . as 
do all of the SZA deep fields. For these fields we use 
NVSS to identify positions for radio sources that are not 
detectable in the SZA 30 GHz data. 

3. ANALYSIS 
3.1. Likelihood Analysis 

As an interferometer, the SZA directly measures the 
Fourier components of the sky brightness. Since it is 
precisely the variance of these components that the power 
spectrum describes, interferometers are ideally suited for 
measuring the power spect rum, without the intervening 
stage of map-making (e.g.. lWhite et al.lll999f) . 

We have used a maximum likelihood method in order 
to extract meas urements of t he po wer spectrum from our 
data. Following lBond et al.l (|1998l ) , we construct a vector 
A of visibilities and a likelihood estimator that is a func- 
tion of A and the CMB power k, assumed to be constant 
across the range of angular scales probed by the SZA. 
We expect our data, and thus the likelihood function, to 
be Gaussian, so that the likelihood is given by: 



C(C) 



1 



• exp 



1 



A T C- X A 



fl 



(2tt) n / 2 \C\ 1 / 2 

where N is the length of the data vector A, and the 
covariance C of the data is written as the sum of a con- 
tribution from the CMB, the diagonal instrument noise, 
and a set of constraints (see below): 

C = K C CMB + Cnoise + C 'constraints ■ (2) 

If our visibilities were independent of one another, the 
contribution to the covariance matrix from the CMB sig- 
nal, kCcmb , would be the identity matrix times the level 
of CMB power, k. We assume a fiat band power, i.e., 
constant £(£+ l)Cg/2w. However, in practice, visibilities 
can be correlated, depending on the separation of the 
Fourier-space (u-v) coordinates they sample. A visibility 
corresponding to a given u-v coordinate is in fact an in- 
tegral over a small patch of the Fourier plane, weighted 
by the antenna's aperture autocorrelation function; the 
overlap integral of this function centered on two neigh- 
boring u-v coordinates yields the correlation between ad- 
jacent measurements. We calculate these overlap inte- 
grals for each of our visibilities, and insert the resulting 
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matrix into the total covariance matrix as indicated in 
equation ((2]). Visibilities that are more than 95% corre- 
lated are binned together. 

The instrumental contribution to the noise, C no i se , 
dominates the weak signal from the CMB anisotropy 
and must therefore be accurately determined. We es- 
timate the noise variance from the scatter of visibilities 
taken on 20 second timescales within a 4 minute scan. 
From the 10 visibilities within each scan, we subtract 
the mean of the real and imaginary parts from each sam- 
ple and then compute a single sample variance from the 
20 mean-subtracted numbers. All averages within the 
analysis are uniformly weighted, as weighting by vari- 
ances derived from so few samples introduces substantial 
biases. The resulting noise estimates are independent of 
one another, so C no i se in equation © is diagonal. We 
have verified through simulation that variances measured 
in this way produce an unbiased estimate of the noise, 
and, hence, of the CMB power level. 

3.2. Constraining Radio Sources 

We expect the sky signal to be dominated by radio 
sources that are pointlike at the SZA's arcminute res- 
olution, principally active galactic nuclei emitting syn- 
chrotron radiation at r adio frequencies. Through the use 
of constraint matrices (jBond et al.lll998HHalverson et all 
2002), we can eliminate the modes of the data that are 
corrupted by these sources, provided that we know their 
positions. The constraint matrix is computed by writing 
a template vector of visibilities that describes the contri- 
bution of the contaminant in question to each u-v compo- 
nent; in the case of a point source at position (I, m) with 
respect to the interferometer's phase center, the contri- 
bution Vp t to the visibility measured with component 
(uj , Vj ) is given by: 

V*t = Ae 2nt(u > l+V > m) , (3) 

where A is an unknown amplitude. By forming the 
outer product of this vector with itself, C l ^ t = V^ t VL, 
we find the covariance created among our N visibilities 
by a source at this position; the amplitude of this co- 
variance remains an unknown, but its shape is fixed by 
the (known) source location. Each constraint matrix is 
then multiplied by a large prefactor (as large as possible 
without causing the matrix to become poorly conditioned 
for inversion) and added to the total covariance matrix 
in our likelihood function, effectively setting to zero the 
weight of the mode corresponding to the source. The un- 
known amplitude is therefore unimportant; the operation 
is equivalent to marginalizing over the source's unknown 
flux. We form such a matrix for each source that we 
wish to eliminate from our data; because each matrix is 
created from a single vector, each has rank one and elim- 
inates one degree of freedom from the N that we have 
measured. Radio sources in our data may have a wide 
range of spectral indices and cannot in general be ex- 
cluded by a single constraint. Instead, we project three 
modes for each source: one flat spectrum and two of dif- 
ferent spectral indices. We therefore lose three degrees of 
freedom for each source projected out. A linear combina- 
tion of these components can approximate any spectral 
index, and the coefficients in this linear combination are 
irrelevant for our technique. 



We remove a component common to the fields within a 
group of four (ground or antenna crosstalk) using a sim- 
ilar constraint. This constraint identifies each visibility 
among a group of four with its three matches in differ- 
ent fields, all of which are measured at the same antenna 
position. Projecting this constraint is equivalent to sub- 
tracting the mean of the four visibilities from each of 
them, and therefore reduces the sensitivity by 25%. 

The radio source positions are derived from several 
data sets. The brightest sources at 30 GHz are iden- 
tifiable directly in the SZA data. For purposes of source 
detection, we combine the long and short-baseline data 
(typical noise of ~ 150 /uJy), providing 30% better sensi- 
tivity to compact sources than would be obtained using 
only the short baselines from which our anisotropy mea- 
surement is determined. A total of 72 sources above 4a 
were detected at 30 GHz, of which ~ 85% have counter- 
parts in the 8 GHz VLA data. Above a significance of 
5<t in the 30 GHz data, 45 of 47 detected sources have 
counterparts at 8 GHz; the remaining two sources were 
found to be extended even at the 20" resolution of the 
long SZA baselines, and are likely heavily resolved by the 
VLA in its C array configuration (3" resolution) . 

As sources below our 30 GHz detection threshold can 
still contribute significant anisotropy power, we also con- 
strain sources detected in the 8 GHz survey. We ex- 
amined several combinations of detection significance 
thresholds in the 8 GHz and 30 GHz catalogs, and con- 
clude that neither catalog is adequate on its own. As 
can be seen in Tabled the measured power continues to 
decrease as the detection threshold in either catalog is 
lowered, while projecting sources from both catalogs re- 
duces the power below that which can be achieved with 
either catalog alone. Additionally, the 30 GHz catalog 
contains sources with no counterparts at 8 GHz, imply- 
ing a population of sources at 30 GHz too inverted, or 
variable, to be detected in the 8 GHz survey. We find 
good agreement in the power resulting from the 5a /6a 
and 4a /5a (30/8 GHz) thresholds and we adopt the lat- 
ter for our analysis. The two fields without 8 GHz cover- 
age are excluded from the analysis shown in Tabled but 
are included in our final result; for these we constrain 
all NVSS sources. We have verified that excluding these 
fields does not significantly change our result. 

For each field, we assemble the noise and constraint 
components of the covariance matrix and vary the level 
of CMB power, evaluating the likelihood of the data ac- 
cording to equation |T]) over a range of values for the 
power k. Treating our fields as independent samples, 
we take the product of each of their likelihood curves to 
form a global likelihood for the experiment; we report the 
maximum of this curve as the most likely power, and the 
points which enclose 34% of the curve's area on cither 
side of the maximum define the 68% confidence interval. 
The likelihood is allowed to extend below zero power. 

We tested this analysis pipeline extensively with simu- 
lated data, including tests in which the data generation 
and power spectrum analysis were performed by different 
parties, using independent software packages, and found 
consistent results in all cases. 

3.3. Jackknife Tests 

Processing the data as described above, we can com- 
pute the relative likelihood of different levels of CMB 
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TABLE 2 

Power Measured Using Various Radio Source 
Catalogs 



30 GHz SZA 


VLA 8 GHz 


Power 


Number of Sources 


Threshold 


Threshold 


[uK 2 ] 




6ct 
5a 
4a 
4cr a 




„ _, -4-79 

1301™ 


33 
41 
64 
53 




20ct 
15a 
12a 


315±^ 
2361S* 
215l|| 


81 
98 
102 




8a 
5a 


110 -54 

105+^ 


111 
180 


5a 
4a 


6(7 

5a 


„+b4 
oo -60 

601S 


171 
244 


Note. — Measurements exclude one 


set of four SZA fields 



that has incomplete VLA coverage at 8 GHz. 

a Only SZA sources with 8 GHz counterparts are constrained 

power, given our data. We now subject the data to three 
jackknife tests to rule out contaminating power from non- 
astronomical sources or from inaccuracies in our noise 
variance estimation. In the jackknife tests, the data are 
split into halves, each of which measures very nearly the 
same set of spatial Fourier components. The data are 
taken in such a way that three of these splits are possi- 
ble: i) a frequency jackknife, wherein we separate even 
and odd frequency bands, as neighboring bands sample 
very nearly the same u-v points, ii) a jackknife in which 
the data are separated into halves by time, and iii) a 
jackknife in which the data are separated into even and 
odd days. The latter two jackknife tests take advantage 
of the fact that the SZA observed the same fields in the 
same way for many days. 

After dividing the data into halves, we difference the 
matching pairs of visibilities. Emission from astronomi- 
cal sources, which is common to both halves, is removed, 
while any contamination that varies with time or fre- 
quency will remain in the differenced (jackknifed) data. 
Note that the requirement that there be matching, un- 
flagged visibilities in both halves for each test results 
in some additional loss of data and degradation of the 
measurement sensitivity, particularly in the two time- 
jackknife tests. We compute variances from the differ- 
enced visibilities, and compute the power in each jack- 
knifed dataset just as for our unjackknifed data, with the 
expectation that the measured power will be consistent 
with zero. The jackknife tests also verify our noise model; 
if the noise in our data has been mis- estimated, we expect 
that the purportedly signal-free jackknifed data will be 
poorly described by our computed covariance matrices. 

The jackknife tests provide evidence that we detect 
a signal correlated with the orientation of our antennas 
relative to the ground. This signal may originate from 
antenna crosstalk or ground emission. As discussed in 
§ 12.11 the observations were designed to allow removal 
of such "ground contamination" by subtracting the aver- 
age of each group of four observations via an additional 
constraint. We find that without removing this compo- 
nent (but including strict cuts on data quality discussed 
in § 12. ip the frequency jackknife test fails badly, with 
a probability to exceed the jackknifed \ 2 of 10~ 5 . After 



TABLE 3 
Anisotropy Measurements for Raw 
and Jackknifed Data 



Jackknife Test 


Power [ijK 2 


PTE 


Frequency 


65l 


-70 

40 


0.19 


First-Second Half 


nt + 106 
zo -81 


0.59 


Even-Odd Days 


-40 


+96 
-49 


0.41 


Unjackknifed Data 


60+ 


65 
Bfi 


0.12 


Note. — Results 


for 4a 


30 


GHz and 



5a 8 GHz source removal and ground sub- 
traction. Two fields without 8 GHz VLA 
data are included using NVSS source con- 
straints. 

implementing the ground subtraction, the three jackknife 
tests all pass, as seen in Table [31 with probabilities to ex- 
ceed the x 2 > 0.1. We therefore employ ground subtrac- 
tion to determine anisotropy power in all the analyses 
presented in this paper. Our analysis procedure is blind 
to the final result, as we have refined our data cuts and 
analysis treatment based only on the jackknife tests. 

3.4. Results 

Repeating the analysis on the complete, ground- 
subtracted, unjackknifed data gives 60l°^K 2 . In the 
next section we estimate the various contributions to this 
power. 

4. ANISOTROPY CONTRIBUTIONS 

To set constraints on secondary CMB anisotropy, we 
must estimate the residual power contributions to the 
SZA measurement from other astronomical sources. In 
this section, we first consider contributions from primary 
CMB anisotropy, undetected radio sources, and diffuse 
Galactic emission, and we end with our constraint on 
the level of secondary anisotropy. 

We note that contributions from undetected radio 
sources and diffuse Galactic emission are in principle dis- 
tinguishable from the CMB by their spectral signature 
across the 25% fractional bandwidth of the SZA data. 
However, it is not possible to distinguish a preferred spec- 
tral index given the low signal to noise of the data. 

4.1. Primary CMB 

At a multipole of I = 4000, the WMAP5 cosmology 
(|Komatsu et al.l I2008D predicts no significant primary 
CMB anisotropy (Figure [T]). However, our window func- 
tion has limited sensitivity to the larger scales at which 
the primary CMB signal is strong. We estimate this con- 
tribution and its variance with simulated observations. 
We generate CMB skies according t o WMAP5 cosmolog- 
ical parameters using CMBFAST dSeliak fc Zaldarriaeal 
1996). We use 35 sets of 43 fields with identical noise 
realizations but different CMB realizations and sample 
them according to the SZA u-v coverage and noise. From 
these simulations, we determine that the primary CMB 
contributes a mean power of 26 ± 5 /iK 2 to our measure- 
ment. Integrating the CMB primary power spectrum 
multiplied by our window function yields a very similar 
estimate. 
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4.2. Undetected Radio Sources 

To assess the level of residual fluctuations from sources 
that are undetected at 8 GHz and 30 GHz at our 5er/4cT 
thresholds (see § El, we simulate 30 GHz source pop- 
ulations using two independent estimates of the source 
density as a function of flux. The first of these is de- 
rived from the SZA blind cluster survey (Muchovej et al. 
2009, in prep.), which provides source counts down to 
1 mjy; we extrapolate to lower flux levels according to 
a power law. This survey is the deepest available at this 
frequency. We also use the source density estimated by 
Ide Zotti et al.l (|2005| ) from an analysis of lower-frequency 
surveys. To assign 8 GHz fluxes to our simulated popu- 
lation of sources we use the spectral index distribution of 
Muchovej et al. (2009, in prep.) measured for 200 sources 
between 5 and 30 GHz. We simulate SZA observations 
with these source populations, applying the source rejec- 
tion procedure described above. We find that residual 
power from undetected sources is small, changing the 
power measured in simulated fields by 20 ± 28 /7K 2 . 

4.3. Galactic Synchrotron Radiation and Free-Free 
Emission 

At lower frequencies, emission from the Galaxy via syn- 
chrotron and free-free radiation contribute significantly 
to the sky brightness , as can be seen in the WMAP 
22 GHz all-sky map (jGold et al.ll200l . We lack sen- 
sitive, high-resolution measurements of these compo- 
nents and are forced to estimate their contribution from 
lower- resolution maps and forecasts derived from them. 
iTegmark et al.l (|2000l ) used degree-resolution templates 
to predict that this emission should contribute less than 

1 /iK 2 at I ~ 4000. Renormalizing their model us- 
ing an analysis of for egrounds measured with WMAP 
(jTegmark et al.l l2004| ). we estimate a contribution of 

2 /iK 2 at 30 GHz from synchrotron and free-free emission 
combined. 

4.4. Galactic Spinning Dust 

A final source of potential foreground contamination in 
our power spectrum measurem ent comes from Galactic 
dust. While the predictions of Finkbcin er et al.l (|1999f ) 
for thermal dust emission at 30 GHz imply < 1 /iK 2 
contribution, there has been evidence that dipole radi- 
ation from small, spinning dust grains may dominate 
the emission at 30 GHz. Such emission is expected to 
peak at frequencies of tens of GHz and to be tightl y 
correlated with Galactic dust (jDraine fc Lazarian]|1998f ). 
Dust-correlated emission within this freque ncy range has 
been observed on various angular scales (jLeitch et al.l 
1997t iKogut et a I Il996t Ide Oliveira-Costa et alJ 120021: 
Finkbeiner et al.l l2002t iFinkbeineij l2004t IWatson et all 
2005HDickinson et al.ll2007ft . During the SZA blind clus- 



TABLE 4 
Anisotropy Contributions 



ter survey, increased map noise was indeed observed in 
fields with more obvious IRAS structure. 

To estimate the contribution of spinning dust to the 
SZA measurement, we take the power spectrum of the 
IRAS 100-/1TO maps of our fields, fit it to log 10 Power = 
Al + B, and extrapolate it to the smaller angular- 
scales that we are sensitive to, but to which IRAS is 
not. We then take the highest measured emissivitics 
of dust-associat e d emi ssion at 30 GHz, as tabulated by 
iDickinson et al.l (|2007f ). and scale this power spectrum 



Source 


Power Contribution [/iK 2 ] 


Primary CMB 


26 ±5 


Undetected Radio Sources 


20 ± 28 


Galactic Synchrotron and Free-Free 


2 


Galactic Spinning Dust 


< 16 


Secondary CMB (SZE) 




and upper limit at 95% C.L. 


149 



to match the observations. Multiplying by the window 
function shown in Figure [1] and integrating over £, we 
get a contribution for each of the observed emissivities 
lying in the range fro m < 1 /zK 2 to 16 j iK 2 . Note that 
the highest emissivity (jLeitch et al.lll997t ) is substantially 
higher than other measurements, and so the correspond- 
ing 16 /zK 2 is likely a substantial overestimate of the 
contribution of this foreground to our measurement. 

4.5. Secondary CMB anisotropy 

We now solve for the constraint on secondary CMB 
anisotropy. We account for the expected contributions 
from primary CMB anisotropy and undetected point 
sources. We do not consider the contributions from 
sources of diffuse Galactic emission, since even the con- 
servative estimates discussed above are subdominant to 
the contributions from primary CMB and undetected ra- 
dio sources. 

Subtracting the residual power for undetected point 
sources and including our 20% calibration uncertainty, 
we are left with 40lg2 A*K 2 f° r the combined primary and 
secondary anisotropy result. 

The constraint on the level of secondary CMB 
anisotropy after subtracting the expected contribution 
from primary CMB anisotropy is 14^g3, /iK 2 . Applying a 
prior of positive power and assuming a Gaussian distri- 
bution results in an upper limit to the level of secondary 
CMB power of 149 /iK 2 at 95% C.L. (78 /iK 2 at 68%) . 

The estimated anisotropy contributions and the re- 
sulting constraint on the level of secondary anisotropy 
are tabulated in Table [4j In Figure [1] we plot the 
SZA anisotropy measurement corrected for residual 
power from radio sources, along with previous measure- 
ments of small-scale CMB power from C BI and BIMA 
(|Readhead et al.ll2004l; [Dawson et aT]|2006l respectively). 

5. DISCUSSION 

5.1. Constraints on eg from cluster simulations 

To understand the implications of our measurement 
for the value of a&, we compare to simulated SZE maps. 
Several groups have carried out large-scale cosmological 
structure simulations and converted three-dimensional, 
simulated universes into two-dimensional projections of 
the Compton y parameter, the frequency-independent 
measure of the magnitude of the SZE. We take these 
maps as inputs to mock SZA observations, and attempt 
to estimate the mean level of power that would be mea- 
sured by the SZA, as well as the scatter in these mea- 
surements. 

We use 60 maps, eac h 0.5° x 0.5° , from the hydro- 
dynamical simulation of IWhite et al.l (|2002l ). generated 
with us = 0-9 i and including cooling and feedback. 
iSchulz fc Wh ite (2003) provide 360 maps from dark mat- 
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ter simulations with erg = 1.0 and gas paste d into clus- 
ter ha loes after the simulation is complete; IShaw et all 
(2008) follow a sim ilar method, produc ing 100 maps with 
cr 8 = 0.77. F inally. IHolder et all (120071) use the Pinocchio 
algorithm of iMonaco et al.l (|2002l ) to generate halo dis- 
tributions and merging histories, resulting in 900 maps 
over a range of er 8 . 

For each set of simulated maps, we form 50 groups of 43 
fields; we pick the maps from the available set randomly 
with replacement. For each map, we simulate SZA ob- 
servations, reproducing the u-v coverage and noise prop- 
erties of the actual data. We process these mock ob- 
servations just as we do the real data, and produce a 
maximum likelihood power for each set of 43 fields. We 
compute the mean and scatter of these powers from the 
50 independent realizations. 

Results for each set of simulated maps are shown in 
Figure [2j The strong dependence of the detected power 
on the value of as is clear, as are significant systematic 
differences among simulations. The maximum likelihood 
power for the secondary CMB anisotropy measured by 
the SZA is shown by shaded horizontal bar in the fig- 
ure. The width of the bar indicates the ±ltr uncertainty 
due to the absolute calibration and the subtraction of 
power from primary CMB anisotropy and undetected ra- 
dio sources. 

Figure O shows that the range of as values consistent 
with the SZA anisotropy measurement is large, reflecting 
significant differences among the simulations. The figure 
shows there is no tension between the SZA anisotropy 
measurement with values of as in the range 0.7 — 0.8 
determined from recent measurements using other tech- 
niques, while values higher than 0.9 appear highly un- 
likely. 

We note that the power spectra derived from simu- 
lated y-maps have greater sample variance than maps 
of Gaussian noise with the same rms power. While we 
can measure the variance in simulation by measuring the 
scatter of many realizations of our experiment, in the ac- 
tual data we assume Gaussianity in calculating our error 
bars. We are therefore underestimating the sample vari- 
ance in our measurement. Using the simulated y-maps, 
we find that the magnitude of this underestimate has a 
weak dependence on the value of as, but for as ~ 0.8 the 
confidence region of our measurement is likely 1.5 times 
broader than the uncertainties on the Gaussian power 
spectrum measurement suggest. 

5.2. Correlations Between Radio Sources and Clusters 

We investigate the possibility that we may bias our 
measurement low by projecting out radio sources associ- 
ated with the very clusters that provide the SZE signal. 
Because the short baselines have a similar response to 
compact radio sources and clusters, in projecting away 
the contribution of the source we may remove consider- 
able flux from the cluster if only the short-baseline data 
are used. This introduces no bias if sources and clusters 
are randomly distributed with respect to each other, but 
we expect radio sources to be spatially correlated with 
clusters, as i ndeed there is some o bservational e viden ce to 
suggest fe.g-. ICoble eTaLl (12001 . ILin fc Mohrl (l200l and 
references therein) . We find that in a pessimistic simula- 
tion, in which every cluster contains at least one source 
that is projected out using only the short baselines, we 
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Fig. 2. — The mean and scatter of the maximum likelihood 
power of several sets of 50 simulated realizations of the SZA mea- 
surement. Each set uses simulated y-maps made with a specific 
technique and a specified rrg. Four different groups input j/-maps 
were used: IHolder et al.l l)2007T ) are shown for multiple values of ag 
as soli d dots, while t riangles indic ate the simulations oflShaw et al.l 
ll2(3u8l) . IWhite et al.l J20OT) and lSchulz fc Whitel (12001) with cr 8 = 
0.77, 0.9, and 1.0, respectively. The points at its = 0.9 and 1.0 are 
offset slightly from those values for clarity. The error bars reflect 
the scatter of the realizations, and thus include the additional sam- 
ple variance from the non-Gaussianity of the y-maps. The dashed 
line and shaded region represent the maximum likelihood power for 
the secondary CMB anisotropy measured by the SZA, including the 
±lcr uncertainty due to the absolute calibration and the subtrac- 
tion of power from primary CMB anisotropy and undetected radio 
sources. 

recover only one-quarter of the actual SZE power in the 
fields. 

Although the short baselines are unable to discriminate 
between radio sources and clusters, the long baseline- 
data are well-suited to making this distinction. On spa- 
tial scales probed by the long baselines, signal from 
a cluster is exponentially damped with respect to the 
short-baseline response, while signal from a point source 
remains constant. In simulations, we confirm that where 
source-cluster correlations bias our measurement low, in- 
cluding the long baselines in the analysis reduces the bias 
by a factor of three. 

In practice, we include the long baseline data in the 
analysis by adding a second bin containing zero cluster 
signal at I ~ 30000, and tying the two bins together with 
constraints for the radio sources. As this procedure in- 
creases the size of the covariance matrix beyond what can 
be manipulated in a reasonable amount of time, adding 
this bin requires further compression of the data by aver- 
aging adjacent scans. Much of the averaged data remains 
> 90% correlated, and we discard any data that is not 
(~ 15%). We have checked that this additional compres- 
sion does not bias the power detected in a single bin in 
^-space. The measured power when including the long- 
baseline data is 50I44 [iK 2 . The consistency between the 
power measurements with and without the long-baseline 
data suggests that source-cluster correlations do not sig- 
nificantly decrease the measured power. 

5.3. Comparison with CBI and BIMA 

Previous measurements at 30 GHz have suggested a 
high-£ excess over the primary CMB anisotropy at a level 
in apparent conflict with the results of this paper, and 
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implying values of erg inconsistent with other contempo- 
rary measurements of the parameter. To explore whether 
the higher excess power could instead be consistent with 
residual power from radio sources, we examine two addi- 
tional prescriptions for projecting sources from the SZA 
data. We also note that accounting for the non-Gaussian 
sample variance discussed above reduces the significance 
o f the discrepancies betw een the measurements. 

iReadhead et~aTl (|2004fl used NVSS to identify sources 
for removal from the CBI 30 GHz anisotropy measure- 
ment, with a limiting 1.4 GHz flux of 3.4 mJy. They es- 
timate the residual contribution from undetected sources 
in the CBI data at ~ 20% of their highest-^ bandpower, 
or ~ 90 /iK 2 . A more precise estimat e of the CBI residual 
power is given 

by[M 

ason et ail (|2003f k for their high-£ bin 
they estimate 115 /iK 2 ±57 /iK 2 . We repeat the SZA anal- 
ysis using the CBI prescription for identifying sources to 
constrain, i.e., sources with NVSS fluxes greater than 
3.4 mJy and with no reference to our 30 GHz source cat- 
alog. This test resulted in a measurement of 378 ±87 /iK 2 
after correction for primary CMB, an excess of 364 /iK 2 
over our nominal result. If we attribute this excess to 
residual point source power and extrapolate it to the £ e jj 
of the CBI measurement according to the respective win- 
dow functions of the two experiments, we find our esti- 
ma tes for the CBI resid ual power to be 189 ± 45 /iK 2 for 
the iMason et all (12003) result and 114 ± 27 /iK 2 for the 
IReadhead et alj (|2004D result. These indicate roughly 
50% higher residual source contamination than estimated 
b y CBI. 

iDawson et al] (|2006h used a 4.8 GHz VLA survey, 
with depth comparable to our 8 GHz survey, to gener- 
ate source constraints for the BIMA experiment. They 
estimate negligible power for residual radio sources in 
their measurement. We approximate applying their se- 
lection criteria to the SZA data by constraining only 
sources detected at 5<r in our 8 GHz data, again ignor- 
ing our 30 GHz catalog. We find a significant increase 
in power, to 107 ± 65 /iK 2 . After subtracting the pri- 
mary CMB contribution this implies 81 ± 65 /iK 2 sec- 
ondary anisotropy, or 67 /iK 2 higher than we measure. If 
wc attribute this instead to residual point source power 
and scale to the BIMA effective I, we would predict 
96 ± 93 /iK 2 of residual radio source power for the BIMA 
measurement. Subtracting this residual power from their 
measurement gives 124^152 A*K 2 secondary CMB power, 
consistent with our result. 

The analysis above suggests that both CBI and BIMA 
may have underestimated the residual contribution from 
undetected radio sources. By contrast with CBI, the 
8 GHz VLA survey of the SZA fields provides a sensi- 
tive measurement of the source population at a frequency 
close to the observing frequency, permitting removal of 
30 GHz power from a typical radio source with a ~ 0.7 
(S v oc v~ a ) down to ~ 15 /iJy, about one tenth of the 
noise level in the SZA CMB fields, and from the most 
inverted sources down to ~ 0.6 mJy, comparable to our 
detection threshold at 30 GHz, and roughly an order of 
magnitude lower than the 30 GHz detection threshold for 
CBI. 

Several of the SZA sources detected at 30 GHz lacked 
counterparts at 8 GHz, implying that they have spectral 
indices more inverted than those we measured between 
5 and 30 GHz (a < — 1 in the most extreme case), or 



were not detected at 8 GHz when the VLA survey was 
obtained due to variability. The former is consistent with 
the expected bias toward more inverted spectral indices 
with decreasing flux density, while the latter is consistent 
with the properties of inverted-spectru m sources, which 
are k nown to be highly variable (e.g., iToffolatt i et alJ 
2005). As can be seen in Tabled! excluding these 30 GHz 
sources from the analysis leads to higher detected power, 
indicating inverted sources near the SZA 30 GHz detec- 
tion threshold are significant contributors to the residual 
source power. They are also the mos t difficult population 
to identify. Indeed, IToffolatti et al.l (|2005j ) demonstrate 
that with slightly different assumptions about this popu- 
lation derived at lower flux levels, a majority of the excess 
power detected by CBI could be explained by residual 
source power. 

6. CONCLUSIONS 

We present results from 30 GHz measurements of 
the CMB angular power spectrum with the Sunyaev- 
Zel'dovich Array (SZA), on scales where the secondary 
anisotropy from the Sunyaev-Zel'dovich effect (SZE) is 
expected to dominate. For a broad bin centered at mul- 
tipole 4066, we find 60 /iK 2 , of which 26 ± 5 /iK 2 is 
the expected contribution from primary CMB anisotropy 
and 20 ± 28 /iK 2 is the expected contribution from un- 
detected radio sources. The resulting constraint on sec- 
ondary anisotropy is 14lg3, /iK 2 , implying an upper limit 
of 149 /iK 2 at 95% C.L. (78 /iK 2 at 68%). 

The SZA results indicate lower secondary anisotropy 
power than previously reported by experiments at 
30 GHz. We show that this discrepancy can be par- 
tially attributed to different methods of identifying the 
contaminating radio sources. The SZA observations em- 
ployed a hybrid configuration of antennas that includes 
both short and long antenna separations, providing high 
sensitivity to arcminute anisotropy while simultaneously 
enabling the detection of radio sources at much higher 
resolution. Accounting for radio sources detected in the 
30 GHz SZA data, and those found with deep VLA ob- 
servations of the SZA fields at 8 GHz, results in the low 
anisotropy power measurement. 

The hybrid configuration was also used to check 
whether anisotropy power was suppressed by removal of 
radio sources correlated with clusters. No significant ef- 
fect was found. 

We show that the level of SZE anisotropy power im- 
plied by the SZA measurement is in good agreement 
with expectations based on simulated y-maps made with 
og ~ 0.8, but in conflict with values of erg of 0.9 and 
higher. The differences among various simulations, how- 
ever, prevent a more quantitative determination of erg. 
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